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HIGH SPEED TRANSPORT CRUISE DRAG 


Abstract 

This report provides scaling laws for the cruise aerodynamics of high speed transport 
wings based on the results of Navier-Stokes computations. Expressions for the various drag 
components are found, together with the corresponding value of (£) m for various values of 
the geometric parameter f which allow for simple optimization of the wing configurations 
with respect to span. It is found that linear theory expressions can be used for this purpose 
provided the coefficients of these expressions for C D and (£)„ available using Navier 

Stokes results. 
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p 


mean value 


t = volume coefficient (= .039) 

Introduction 

An improved understanding of the cruise aerodynamics of the high speed transport 
wing requires that both theoretical and experimental methods be applied to define al- 
ternative configurations suitable for the high speed flight regime. Previous aerodynamic 
analysis has relied heavily on linear analysis together with extensive testing of cruise con- 
figurations. In order to prepare for future designs new methods of analysis utilizing modern 
CFD techniques appropriately validated by wind tunnel experiments must be applied be- 
fore configurations are selected for experimental development in large facilities. 

The approach used for this report utilizes the results of Navier Stokes computations 
for variations primarily in j about a Boeing wing used as a baseline. The lift and drag 
increments are assumed to have the same dependence on the planform parameters as 
suggested by linear theory but the dimensionless coefficients are now evaluated using the 
non-linear Navier-Stokes results. In this way a simple general expression for L/D for the 
configurations can be formulated. 

Computations have been made to demonstrate the approach and results are given for 
(L/D) m permitting a reduction in the experimental testing required. An upper limit to 
(L/D) m is found and the corresponding optimum value of f is given. 

Analysis 

Linear theory suggests that the drag coefficient of a supersonic configuration can be 
written in the form 


C D = C DF + fr V (^ Ka + ±Cl^ 


K w + 2 0 2 K, 



( 1 ) 
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where 


Cdf is associated with skin friction drag 
K 0 is associated with wave drag due to thickness 
K v is associated with vortex drag 
K w is associated with wave drag due to lift 
Cl is the lift coefficient 
and 

P = Ml- 1 


In the Navier-Stokes computations f was varied by a scale factor f from a nominal of 
f = .4336895, ie f = .4336895 f, 1.25 > f > .625, and p was varied from a nominal value 

of p = .2332929 to a value of p = .2210521 

A general equation for the drag coefficient was postulated of the form 

C D = A + Bp 2 / 2 + 3E(C + D p) (2) 

where A,B,C and D are coefficients to be evaluated using the Navier Stokes results and 
a mean value of p given as p = .2281858, is used. The expression (2) is analogous to (1) 
and is useful only if the value of C D that results is sufficiently close to that given by the 
Navier-Stokes results. 

(1) Zero Lift 

First mean value of the drag at zero lift were found from the Navier-Stokes computa- 
tions for each value of f: 
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Cdo = -008546 (mean of wings 15, 25, 35 and 42;/ = 1.25) 

= .007541 (mean of wings 1, 2, 3 and 4; / = 1) 

= .006700 (mean of wings 16, 21, 31, and +3;/ = .875) 

= .005832 (mean of wings 12, 22, 32, and 44;/ = .75) 

= 005176 (mean of wings 13, 23, 33, and 45/ = .625) 

and the constant A and B evaluated by matching the results for f = 1 and f = .625; this 
gives 


A = .003660 and B = .074537 


Thus 


C Do = .003660 + .074537p 2 / 2 (3) 

Table 1 shows the values of Cdo computed, for each wing, by the Navier Stokes equa- 
tions with the result given by equation (3). 

Table 1 - Cdo 


Navier-Stokes 


Mean 

Equation 3 

Wing 15 

.008854 

.008846 

.008572 

Wing 25 

.008551 



Wing 35 

.008636 



Wing 42 

.008342 



Wing 1 

.007552 

.007541 

.007541 

Wing 2 

.007734 



Wing 3 

.007478 
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Wing 4 

.007402 



Wing 16 

.006851 

.006670 

.0066631 

Wing 21 

.006144 



Wing 31 

.006608 



Wing 43 

.006577 



Wing 12 

.005799 

.005532 

.005843 

Wing 22 

.005845 



Wing 32 

.005851 



Wing 44 

.005833 



Wing 13 

.005151 

.005176 

.005176 

Wing 23 

.005168 



Wing 33 

.005180 



Wing 45 

.005205 



(2) Non Zero Lift 





Similarly, the mean values of the quantity Pe-fe ^L were found from the Navier Stokes 
computations at the value of C D corresponding to that for the maximum computed value of 
(L/D); these values at / = 1 and / = .625 were used to find CandD (since Cn cf n ° { = C + D f*) 
thus 

Cp - C ool _ 19506Q (mean of wings 15 25 , 35 and 42; / = 1.25) 
p 

= 2.6683 (mean of wings 1, 2, 3 and 4; / = 1.0) 

= 2.34708 (mean of wings 16, 21, 31 and 43; / = .875) 

= 2.08039 (mean of wings 12, 22, 32 and 44; / = .75) 

= 1.81567 (mean of wings 13, 23, 33 and 45; / = .625 
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and the constants C and D evaluated by matching the results for f — 1 and f .625, this 


gives 


C = 1.2691 D = 1.3992 


so that 


£R — £doL = 1.2691 + 1.3992/ 2 (4) 

Cl P 

Table II shows the value of computed for each wing by the Navier-Stokes 

u l P 

equations with the result given by equation (4). 

Table II - 

C L P 


Navier-Stokes 


Mean 

Equation 4 

Wing 15 

2.88637 

2.95060 

3.0399 

Wing 25 

2.96967 



Wing 35 

2.97739 



Wing 42 

2.96598 



Wing 1 

2.74649 

2.66834 

2.6683 

Wing 2 

2.59041 



Wing 3 

2.66748 



Wing 4 

2.66903 



Wing 16 

2.28755 

2.34708 

2.34036 

Wing 21 

2.35686 



Wing 31 

2.37203 



Wing 43 

2.37183 
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Wing 12 

2.08080 

Wing 22 

2.07449 

Wing 32 

2.08409 

Wing 44 

2.09281 

Wing 12 

1.80017 

Wing 23 

1.80800 

Wing 33 

1.82395 

Wing 45 

1.83057 


2.08039 2.05615 


1.81567 1.81567 


Total Drag 

The expression for total drag is therefore 


C D = .003660 + .074537p 2 / 2 + [1.2691 + 1.3992/ 2 ] 


( 5 ) 


Comparing (1) and (2) it can be seen that 




— = 1.5989 
512 


K v = 2t r 



C = 3.4582 


and 

K. = £ (7) D = 21293 

for (3 2 = (2.4 2 — 1) 

Thus good agreement is found between equation (1) and the Navier Stokes results if 
(with t = .039) 


K e = 1 598; K v = 3.4582 
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and 


K w = 2.1293. 

Equation (1) can then be used in the analysis of the aerodynamics of the wing. 

With Cq in the simple form the maximum value of L/£) for given f is found by varying 
Cl and is 


L_ 

D, 


which occurs at a lift coefficient Cl given by 


-I 

\(A + Bp^P)(C + Df)p] 

L 2 

f J 




( 6 ) 


(Cl) ) m = 


(A + Bp 2 f 2 )p 


f(C + DP) 


Figure 1 shows as a function of f for the family of configurations. 

An optimum value of / can be found at which ^ j is a maximum by differentiation 
of 6 with respect to /, thus the optimum occurs at / given by 


2 / 2 / 1 „ 
+ -zr L = 0 


P + a P + c f 
when a = , c = § which can be solved for / Q pt to give 


/opt = [( fl2 + ° 2 + 14ac ) 2 “ ( a + c )] 

= .5552 for p = .2281858 

Thus the optimum value of / is substantially less than 1 (the baseline value) and 
corresponds to a value of j = .5552 x .4336895 = .2408. 
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If i >> .2408 the wave drag due to thickness (which varies as (7) ) and the wave drag 
due to lift (which varies as (f) are too large. If f « -2408 the vortex drag (which varies 
as (f) *) is too large. 

The value of ( with / = .5552 and p = .2281858 is given by equation (6) as 
V ” / m 



This is the maximum value of (£) m that can be attained by this family of wings (and 
compares with (j>) m = 7.405 for / = 1). 

Concluding Remarks 

It is found that the expression for supersonic cruise drag of a wing given by linear 
theory can be used and the components of drag identified by evaluating certain coefficients 
in this expression using the results of Navier-Stokes computation. This then facilitates the 
calculation of the maximum % attainable by this family of wings and the optimum f at 
which this occurs. 

Thus it is found that the supersonic cruise drag can be expressed as 


Co = C DF + 5HrV (y) 2 K + icl-^ [k, + 2 KJ? (f ) : 2 


where K 0 ,Kv and K w are evaluated as K 0 = 1.5987; K m = 3.458 and K w = 2.129. 

The maximum value of (^) m = 8.58 was obtained at a value of f = .2408 suggesting 
that a smaller value of f than the nominal value of .4337 gives improved performance. 
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